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Recent measurements of selected QCD processes at the Tevatron are reviewed and con- 
fronted with theoretical calculations. Results on inclusive jet production at large transverse 
momentum (pr) are compared to predictions from next-to-leading order (NLO) perturba- 
tive QCD (pQCD). Kinematic distributions of jets with light and heavy flavor produced in 
association with electroweak bosons are compared to expectations from leading-order (LO) 
QCD calculations supplemented with parton-shower models. Properties of QCD radiation 
in hard-scatter events are investigated using azimuthal correlations between two leading 
jets in multi-jet events, while aspects of softer radiation are examined through properties 
of energy flow within jets. Finally, characteristics of soft interactions underlying the hard 
scatter are explored in the context of tuning the parameters of phenomenological models 
employed in QCD Monte Carlo event generators. 



The ongoing physics program at the Tevatron offers an unprecedented opportu- 
nity to explore a wide variety of QCD issues at a level of high precision. In Run II, 
the collision energy increased to y/s — 1.96 TeV, from 1.8 TeV in Run I, providing 
an increase by a factor of 5 in the inclusive jet cross section at pt~ 600 GeV/c. 
Now, in the early phase of Run II, the CDF and D0 experiments have already 
collected substantially more luminosity than in all of Run I. As a consequence, the 
current reach in px of the inclusive jet spectrum exceeds that of Run I by ~ 150 
GeV/c, and events with dijet masses in the range of 1.2-1.3 TeV/c 2 have been col- 
lected by both collaborations, probing interaction length scales of 10 -19 m. With 
steadily increasing data samples, we are all looking forward to improved precision 
measurements and to exploiting the discovery potential for new physics. 

Production of jets at large pr has been long recognized as one of the most in- 
teresting probes of pQCD. Within the pQCD framework, hard-scattering processes 
are described by a convolution of partonic cross sections with parton distribution 
functions (PDFs). Thus, measurements of high-pr interactions probe PDFs and 
the strong coupling constant, a Sl and test current phenomenological tools (among 
others, high-order matrix-element calculations, models of parton showers, the in- 
terplay between these two approaches, resummation techniques etc.). Deviations 
from pQCD, not accommodated by uncertainties in these components, may indi- 
cate presence of physics beyond the Standard Model (e.g., parton compositeness). 

Production of objects at low and moderate px is sensitive to fragmentation, 
multiple-parton interactions, and to the nature of remnants of the colliding par- 
ticles not participating in the hard collision. Studies of these effects help develop 
phenomenological descriptions of such soft physics and in the tuning of event gen- 
erators. 

The hard-scattered quarks and gluons in the final state fragment into jets of 
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particles that subsequently interact in the detectors. Experimentally, jets are usu- 
ally defined by their energy deposition in multiple calorimeter towers. For Run II, 
both CDF and D0 collaborations adopted a new iterative fixed-radius cone algo- 
rithm. While Run-I jets consisted of i?T-weighted sums of calorimeter towers, the 
Run II cone algorithm defines jets by summing 4-momenta of towers within a cone, 
using midpoints between such jets as additional seeds for the algorithm to minimize 
sensitivity to infrared divergence of gluon radiation £Q. Analogous procedure can 
be applied at parton or particle levels. (CDF is also using their Run I version of the 
cone algorithm JetClu [2] for easier comparison with previous results.) Both col- 
laborations continue using the Kx clustering algoritm, based on relative transverse 
momenta between clustered objects. 




Pr [GeV/c] Pr [GeV/c] 

Fig. 1. Left: Inclusive jet cross section measured by D0 in three rapidity bins using the 
midpoint algorithm vs. jet pr- The error bars are statistical only. Right: Ratio of the cross 
section for the central rapidity bin to NLO pQCD; the outside band represents the current 
experimental systematics, while the inner grey band shows uncertainty from PDFs. 

The high- a; behavior of PDFs has been scrutinized intensely for the past several 
years. The Run I measurement from D0 3 of the inclusive jet cross section in 
several regions of pseudorapidity has provided powerful new information for the 
global determinations of PDFs by the CTEQ jt| and MRST |S] groups. Similar 
measurements are being pursued in Run II. D0 has recently presented the prelim- 
inary inclusive jet cross section as function of px in 3 rapidity (y) bins (Fig. 
and a dijet cross section for central rapidities \y\ < 0.5) (Fig. |5J, based on data 
corresponding to an integrated luminosity of 143 pb _1 . Given the large systematic 
errors, the data are in good agreement with NLO pQCD jS], for the parameters 
of the calculations noted in the figures. So far, there is no indication of parton 
substructure or presence of narrow resonances decaying to jets. The uncertainties, 
shown as bands on the plots of data/theory ratio, are dominated by the current un- 
derstanding of the jet energy scale, and are expected to be reduced in near future. 
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Fig. 2. 
Right: 
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Left: Inclusive dijet cross section as function of dijet mass measured by D0. 
Ratio to NLO theory; the bands show the experimental and PDF uncertainties. 



CDF has presented a measurement of the inclusive jet cross section for central 
rapidities (0.1 < \y\ < 0.7) using a data sample of 177 pb -1 . It is compared to 
NLO pQCD and to Run I result in Fig. (both Run II and Run I results 
were obtained with their JetClu algorithm). The expected increase in the Run-II 
cross section as function of pr for the larger y/s is clearly visible. While there is 
reasonable agreement between data and theory, it will be interesting to re-examine 
these comparisons when experimental uncertainties become reduced. 
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Fig. 3. Left: Central inclusive jet cross section measured by CDF using the JetClu 
algoritm vs. jet Et- Experimental and PDF uncertainties are also shown. Right: Ratio of 
jet cross sections for Run II to Run I, compared to NLO pQCD. 

For a somewhat smaller data sample of 145 pb -1 , CDF measured the inclusive 
jet cross section over central rapidities using the Kt clustering algorithm |S] and 
jet-size parameter values of D — 0.5,0.7,1.0. The results are presented in Fig. 0] 
for D = 1.0, in terms of ratios to NLO pQCD 0, using QCD scale p^ ax /2, and 
to the prediction from PYTHIA jHj. While the ratio to NLO exhibits an excess at 
lower values of pt, that grows as the value of D increases (not shown), no such 
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effect is seen when comparing to pythia with parameters tuned to CDF data on 
the properties of the "underlying event" (so called Tune A, see end of paper) . This 
behavior may be due to the influence of particles in the underlying event on jet 
properties; the measured cross section is not corrected for this contribution. The 
underlying-event effects are not present in the NLO calculation, while they are 
included in the phenomenogical model implemented in pythia. 
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Fig. 4. Left: Ratio of the central inclusive jet cross section measured by CDF using the 
Kt algoritm (with D — 1.0) to NLO pQCD, as function of jet pr- Right: Similar ratio to 

prediction from pythia (Tune A). 



Production of electroweak bosons in association with several jets has become 
an area of increasing interest for pQCD. The cross section for processes with n 
jets is directly sensitive to a™, and the presence of the heavy boson is expected to 
assure a more reliable perturbative calculation. NLO predictions are available up 
to n — 2, while LO calculations (alpgen jTH] code being primarily used by CDF 
and D0) exist for larger values of n. Frequently, such LO calculations are interfaced 
at the parton level to event generators such as pythia or herwig JT] to include 
the physics of soft emissions, described through the parton-shower approach, and 
to facilitate full simulation of experimental conditions following hadronization and 
modeling of underlying event. Such "enhanced-LO" procedures are of obvious in- 
terest for investigations of many related physics processes, including top-quark and 
Higgs-boson production, and searches for SUSY, for which W/Z + jets production 
represents a major background. Figure shows CDF data on the W+ > n jets 
cross section vs. the inclusive number of jets, and the dijet invariant mass distri- 
bution in W+ > 2 jet events, for jets with Er > f5 GeV and \rj\ < 2.4, for a 
data sample of 127 pb -1 . The measurements are compared to expectations from 
ALPGEN+herwig (including full detector simulation and event reconstruction) for 

two choices of QCD scale of Myy and \J < Pr ^ ct >; the calculations bracket the 
data, but the theoretical uncertainty is large at LO. 

CDF and DO have developed algorithms to identify jets originating from b- 
quarks (6-tagging). Figure illustrates results from DO for 174 pb -1 of data on di- 
jet mass distribution in W+ > 26-tagged jet events, compared to ALPGEN+pythia 
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Fig. 5. Left: CDF cross sections for W + jets production vs. the inclusive number of 
jets, compared to the LO prediction from Alpgen+herwig, for two choices of QCD scales. 
Right: Invariant mass distribution of two leading jets in CDF W+ >2 jet events, compared 
to the same LO calculations. Uncertainty from jet energy scale is shown by the band. 



predictions for several contributing processes, and for jet px distribution in Z + b 
events compared to pythia + experimental background. Good agreement is ob- 
served between data and expectations; forthcoming high- luminosity data will chal- 
lenge NLO calculations with much reduced uncertainties. 
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Fig. 6. Left: Invariant mass distribution of two leading jets in W+ > 2 b-tagged jet events 
measured by D0. Data are compared to the simulated processes using Alpgen+pythia. 
Right: Et spectrum of fo-tagged jets associated with Z production measured by D0. The 
shaded histogram represents background contributions estimated from data. The open 
histogram is the sum of background and simulation of Z + b production using pythia, 

normalized to data. 
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Aspects of radiation in hard-scatter events can be studied without reconstruct- 
ing multiple jets, by simply investigating the azimuthal angle between the two 
leading jets, A^dijet- At lowest-order, the two jets are expected to be produced 
back-to-back, A0dijct~ tt- Thus, the distribution in A0 dijet away from 7r reflects 
the presence of higher-order radiation. D0 data on A</> dijet for 150 pb _1 , normal- 
ized by the dijet inclusive cross sections integrated over the same phase-space, are 
presented in Fig. [7] for four ranges of pr of the leading jet (the second leading jet 
was required to have pr > 40 GeV/c). Data is compared to pQCD LO and NLO 
calculations for 3-jet production using nlojet-I — h \T2\ . It is apparent that the LO 
calculation has a very limited range of applicability, while the NLO prediction, 
with up to 4 partons in the final state, provides a good description over a large 
range of A0 dijet- pythia and herwig employ parton-shower models to describe 
higher-order QCD effects, herwig agrees with the data well over the entire range 
of A</> dijet; the default version of pythia provides a rather poor representation 
(dashed lines). However, the amount of radiation in PYTHiA-generated events can 
be adjusted, e.g., by varying the scale factor for the maximum p? allowed in the 
initial-state parton shower (PARP(67)). The bands correspond to the variations in 
this parameter between the default value of 1.0 and 4.0 (used in CDF Tune A). 
The higher values provide a better description, and the data clearly can be used to 
further tune the pythia model. 




Fig. 7. Left: A<j> dijet distributions from D0 in four regions of leading-jet pr- The solid 
(dashed) lines show the NLO (LO) pQCD predictions. Right: Comparison of Acfi dijet 
distributions to results from herwig (solid line), default pythia (dashed line), and a 
possible variation in the level of initial-state radiation in pythia (band). 
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The distribution of energy flow within the jet cone is sensitive to soft parton 
radiation. Figure 03 shows CDF data on the fraction of jet pr integrated over the 
outer part of the cone (0.3 < r < 0.7 for R = 0.7 size cones) vs. jet pr, using the 
midpoint algorithm for 170 pb _1 . pythia, with the Tune-A set of parameters, as 
well as herwig, properly describe the transition of jet shapes from gluon-dominated 
at low pr to quark-dominated at large pr', the default version of pythia, and 
PYTHIA with multiple-parton interactions switched off, undershoot the data. 




Fig. 8. Left: CDF measurement of "jet-pr fraction" in the outer part of the jet cone, 
0.3 < r < 0.7, vs. jet pr, compared to prediction from pythia (Tune A). Expectations for 
pure quark and gluon-initiated jets are also shown. Right: Same data compared to results 
from herwig, and pythia with default settings and with no multiple-parton interactions. 

In addition to objects emerging from the hard scatter, events contain soft par- 
ticles produced through several mechanisms, often collectively referred to as un- 
derlying event (UE). These include fragments of beam remnants, multiple-parton 
interactions, and initial and final-state radiation. Not all of these contributions are 
calculable within perturbation theory, and must therefore be modeled using Monte 
Carlo programs tuned to data. Since particles from UE influence jet properties, 
especially at low pr, proper modeling of UE is essential to account for such effects. 
To enhance sensitivity to UE, in Run I CDF studied ^3] distributions of charged 
particles with pr > 0.5 GeV/c and \y\ < 1.0 in the region of the azimuthal plane 
tranverse to the leading jet (defined using either charged particles or calorimeter 
information) . Examples of similar Run II studies are presented in Fig. [5] for the 
density of charged particles in the transverse region vs. pr of the leading jet, and 
vs. pr of the charged particles in the transverse region (in bins of pr ot the leading 
jet). The event generators require tuning of phcnomcnological parameters to data 
to describe the details of such distributions. CDF studies established a special tune 
of pythia parameters, Tune A, which provides a significantly improved descrip- 
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tion of UE properties compared to default settings. Inclusion of other processes in 
the tuning is expected in the near future, and should help further constrain the 
parameters and verify the universality of this approach. 

A prediction from PYTHIA Tune A is that, at ^/s = 14 TeV, «12% of all interac- 
tions will result in hard scatters with px > 10 GeV/c — potentially an important 
effect on measurements in the high-luminosity environment of the LHC! 
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Fig. 9. CDF studies of charged particle density in the transverse region relative to the 
leading jet (see text) vs. pr of the jet (left), and vs. pr of the charged particles in the 
transverse region (right), compared to expectations from tuned pythia (Tune A). 

Further improvements in the understanding of aspects of QCD at the Tevatron, 
and development of related physics-analysis and simulation tools, will undoubtedly 
be of great benefit to the physics program at the start-up of the LHC. 

It is a pleasure to thank my colleagues at CDF and D0 for many helpful discussions. 
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